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Abstract. The experimental determination of the lifetime of pionium provides 
a very important test on chiral perturbation theory. This quantity is determined 
in the DIRAC experiment at CERN. In the analysis of this experiment, the 
breakup probabilities of pionium in matter are needed to high accuracy as a 
theoretical input. We study in detail the influence of the target electrons. They 
contribute through screening and incoherent effects. We use Dirac-Hartree-Fock- 
Slater wavefunctions in order to determine the corresponding form factors. We 
find that the inner-shell electrons contribute less than the weakly bound outer 
electrons. Furthermore, we establish a more rigorous estimate for the magnitude 
of the contribution from the transverse current (magnetic terms thus far neglected 
in the calculations). 



PACS numbers: 36.I0.-k, 34.50.-s, f3.40.-f 

1. Introduction 

In experiments such as atom-atom scattering, atom-electron scattering, nuclear 
scattering at high energy, or photo-production of e + e~ pairs on atoms one faces 
the situation of complex systems undergoing transitions between different internal 
states. For an atomic target, excitation affects the electrons individually. Thus 
the "form factor" for the target-inelastic process takes the form of an incoherent 
sum over all electrons, as opposed to the coherent action of the electrons (and the 
nucleus) in the target-elastic case. From this observation it is immediately obvious 
that the target-inelastic cross section is proportional to Z, the number of electrons in 
the target atom, whereas the target-elastic process scales with Z 2 . Since incoherent 
scattering off the excited target electrons increases the cross section from its value 
due to coherent scattering off the atom, the effect is sometimes referred to as "anti- 
screening" [I, 2]. At large momentum transfer the anti-screening correction can be 
accurately approximated by increasing the coherent scattering cross section by the 
factor (I + l/Z), see e.g. [3, 4, 5]. We will demonstrate in this article, however, that 
this simple re-scaling argument is not sufficiently accurate in general. At the same 
time we will show how to obtain far superior results in a Dirac-Hartree-Fock-Slatcr 
approach with manageable numerical effort. 

As an application of the general formalism discussed in this paper, we put our 
focus on the particular situation pertaining to experiment DIRAC. This experiment, 
currently being performed at CERN [6], aims at measuring the lifetime of pionium, 
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i.e., a 7r + 7r~ pair forming a bound state as an exotic, hydrogen-like atom. While 
moving through matter, this system can breakup via electromagnetic interaction, or 
it can annihilate in the strong hadronic decay n + n~ — > 7r°7r° if the pionium is in an 
s-state. The lifetime for this decay is related to the mr scattering length which in 
turn has been calculated in the framework of chiral perturbation theory. Experiment 
DIRAC therefore provides a crucial test for a theoretical low-energy QCD prediction 
of pion-pion scattering. The pionium is formed from pions produced in the collision 
of a high energy proton beam on a (heavy clement) target foil on condition that 
the two charged pions have small relative momentum. The pionium production rate 
can be inferred from the double inclusive production cross section for 7r + 7r~-pairs 
without final state Coulomb interaction. The kincmatical conditions of the experiment 
(target thickness of 100 to 200 fim, compared to a predicted decay length in the lab 
of roughly 15 to 20 /im) imply annihilation of the pionium through strong interaction 
still within the target foil. However, the experiment does not measure the neutral 
pions. Instead, charged pions are actually detected, and pions from electromagnetic 
breakup of "atomic" pairs are distinguished by small relative momentum and small 
angular separation from a large background of accidental "free" pairs with arbitrary 
momentum and no directional correlation. Comparing the observed "ionized" atomic 
pairs with the number of pionium atoms produced for a given target material and 
thickness, one can extract the lifetime for the hadronic decay Thus very accurate 
cross sections for the electromagnetic interaction between pionium and the target 
material are needed as an input in the analysis of the experiment. 

The paper is organized as follows: In section 2 we review the well established 
formalism for one-photon exchange, applied e.g. in electron- hadron scattering. In 
section 3 the cross sections coming from the transverse photons are estimated with 
the help of the long- wavelength limit. The formulas for the evaluation of atomic 
form factors and scattering functions in the framework of Dirac-Hartree-Fock-Slatcr 
theory are derived in section 4, followed by an analysis of simple alternative models 
in section 5. Although the results presented in section 6 have been obtained in the 
context of experiment DIRAC, many of the conclusions drawn in section 7 are also 
valid more generally in atomic scattering. 



2. Formalism 



In [7] we applied the semiclassical approximation to calculate the coherent (target- 
elastic) cross section for pionium-atom scattering, demonstrating that this part can 
indeed be determined with the desired accuracy of 1%. However, contributions not 
yet included in [7] need to be added to achieve an overall accuracy at the 1% level. 
Here we are interested in those processes, where the atom is excited together with the 
pionium. We will treat the problem within the PWBA. For our derivation we will 
follow closely the standard formalism of electron scattering as found, e.g., in [8, 9, 10]. 

We only calculate the lowest order result. The relevant Feynman diagram is 
shown in figure 1. The cross section for this process is given by 

° = 4lj2^ 2MA2Ma J d q W? ' (1) 

where I is the incoming flux, and and are the electromagnetic tensors 

describing the electromagnetic interaction of atom and pionium with the photon. As 
we are not interested in the specific final state of the atom, we can average over the 



Coherent and incoherent atomic scattering: Formalism and application 



Pn 



P 



(7171) 



n 




Figure 1. The lowest order Feynman diagram for the simultaneous excitation 
of projectile (pionium) and target (atom). The atomic momenta are P A and P' A 
before and after the collision, those of the pionium P n and P^. The momentum 
of the exchanged photon is q = P^ — P n = —{P' A — P\)- 

possible initial spin and sum over all final states and directions corresponding to a 
specific final state momentum P'. For the atom we have 

= E <°< Pa\J^\X, Pk) (X, P A |J"|0, Pa) (2tt) 4 S 4 (Pa -q-Pk) (2) 

and similarly for the pionium if its final state is not resolved either. 

The electromagnetic tensor can only be a function of P and P', or equivalently, 
of P and q. Gauge invariance or current conservation restricts the possible tensor 
structure of W^ v even more. It is a well known result that the electromagnetic tensor 
in this case only depends on two scalar functions W\ and W 2 that are functions of q 2 
and Pq alone. The electromagnetic tensor is then given by 

r + tf) WM + (p- - »f) - ?±f) 2M^),a) 

Since the cross section depends on the product of the tensor for the atom and the 
pionium, we calculate this product in terms of Wi and W2: 

WT^Vn - 3Wi, n Wi,A + ^-1 + Wi,nW 2 , A + (-1 + ^ W 2 ,nW 1>A 

/ luA\ 2 

+ w 2-nW 2 , A , (4) 

where 7 is the (relative) Lorentz factor between atom and pionium, A = — Pxq/Mx, 
oj = P]jq/Mu are (minus) the energy of the exchanged photon in the atom and pionium 
rest frame, respectively. Following the argument of [11, 7] we expect the cross section 
to be dominated by the charge operator (termed "scalar interaction" in [7]). Thus we 
neglect at this point all terms containing a W\ . We note also that W\ vanishes in the 
elastic case for spin particles. However, the contribution to the cross section coming 
from W\ will be discussed in section 3 below in the analysis of the transverse part of 
the current operator j. Keeping only the last term and assuming that the pre- factor 
is dominated by 7 (7 is in the range 15-20 in the DIRAC experiment) we get 

W^W^u w j 2 W 2M W 2A . (5) 
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Although this coincides with the naive estimate obtained by superficially identifying 
the leading power of 7 (assumed to be dominant), one should keep in mind that (as 
will be shown below) the range of q 2 starts at about (w/7) 2 or (A/7) 2 , where the 
last three pre-factors are of the same order. The magnitude of 7 should therefore 
not be mistaken as a justification for neglecting the terms containing W\. A critical 
assessment of the relative importance of W\ as compared to W2 is postponed until 
section 3. 

In our application to pionium-atom scattering the masses Ma and Mn will be 
much larger than the momentum transfer q of the photon. Therefore we will neglect 
recoil effects on the atom and pionium. We can then identify A and to as the excitation 
energy of the atom and pionium in their respective rest frames. Therefore qo and q z 
are fixed by the values of u> and A. In the following we will always denote the spatial 
part of the photon momentum in the rest frame of the atom by k and in the rest frame 
of the pionium by s. In the rest frame of the atom we then have 

Qo,A = -A 

_A w (6) 

and in the rest frame of the pionium 
<7o,n = w 

A LU 

q z ,n = s z — - — . 

7/3 13 

Also q 2 is given by 

,2 



(7) 



/? 2 7 



dcodAdq ± — — — . (9) 



Replacing the integration over d 4 q = ^dw dA d 2 q±, we obtain 

4a 2 W 2 , n (q;,g 2 )Ty 2 , A (A, g 2 ) 

From (8) we note that qf may become negative if for example the pionium is initially in 
an excited state and gets de-excited (w < 0) while the target atom gets excited (A > 0) 
from its ground state. For this process the cross section (9) has a formal singularity, 
as pointed out recently in [12]. Formally the divergence arises from integrating over 
the impact parameter all the way to infinity. Under the conditions of the experiment 
DIRAC, however, this formal divergence is negligible compared to the principal value 
of the cross section integral. 

One of the advantages of this derivation is that the Wi are scalar functions; we 
can therefore evaluate them in the respective rest frames, even though the relative 
motion of atom and pionium is relativistic. We now establish a relation between W2 
and the electromagnetic transition currents. As was already done in [7] we assume 
the charge operator to be the dominant contribution. In the atom's rest frame W2 is 
related to the 'OO'-component of the tensor through 



„4 



w 2 , A = ^W»\, (10) 



Coherent and incoherent atomic scattering: Formalism and application 



5 



see, e.g., [13, 14], with k 2 = k 2 = A 2 — q 2 . Here we have added an index "rf" as 
a reminder that the '00 '-component in the rest frame has to be taken. From the 
definition of W^ u we get 

W ™ = idul £ (°\ J °^l)\X,P A ) ( X > Pa\J O (i)\0) (27r)^ 4 (PA - q P A ). (11) 

Rewriting this expression in terms of the (non-relativistic) density operator and again 
neglecting recoil effects we find 

WT = < K9)I X ' E x ) (X,E x \p(q)\0) S(E , A + A - E x ) (12) 
x 

= Y J \ F ™A<l)\ 2 5{E^A + &-E x ). (13) 
x 

(Replacing the energy E by the rest mass M in these weakly bound systems only 
introduces an error of the order a 2 ). Finally we obtain for the cross section: 



a = / dwdAd 2 q ± - 



/? 2 s 4 fc 4 



^I^X0,A(fc)| 2 ^0,A + A-^) 



X 



\ F x>o.n(s)\ 2 S(Eo.n + w - S^O 



A"' 



(14) 



As a verification of this formalism we reproduce the result of [7] for target-elastic 
scattering. In this case A = and k 2 = q 2 , and the only possible final state is X = 
(assuming no degeneracy of the ground state) . We find 



d^dV^l|F o, A (g)| 2 



^2\F X oM\ 2 S(Eo,n+u-E x ) 



x 



(15) 



and -F o,a(<z) is the elastic form factor of the atom. This is identical to the equation 
derived in [7] in the Coulomb gauge (fe • A = 0). 

At this point we should add a few comments regarding the gauge invariancc 
and the factor q 4 /k 4 (and likewise q 4 /s 4 ) in (10): In determining the general tensor 
structure of W^ v in (3) the gauge invariance (or alternatively current conservation) 
was used. The magnitude of the component of j along the direction of q is then fixed 
by p. Only the transverse parts of the current remain independent quantities. Our 
result therefore agrees with the one in Coulomb gauge, as in this gauge the component 
of A along the direction of q vanishes. An alternative approach might start directly 
from W^ v without decomposition into W\ and VF 2 - Assuming that in the rest frame 
of either the atom or the pionium W ao dominates, one can approximate the product 
of the electromagnetic tensors by 

^n«AW (16) 

where the factor j 2 comes from the Lorentz transformation of one tensor into the rest 
frame of the other. Now using (16) instead of (10), the cross section reads 



dw dA d 2 qj_ 



4a 2 J_ 

P 2 q 4 



\ F xo,A(k)\ 2 6(E , A + A - E X ) 



L X 



\ F x'o.n(s)\ 2 5(E .n + u - E x ,) 



X' 



(17) 
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This differs from (14) by a factor q 4 /k 4 ■ q 4 / s 4 . In the elastic case (17) corresponds 
to the result of [7] in the Lorentz gauge. However, the approximation (16) and 
therefore also (17) is not gauge invariant, whereas (14) is by construction. We will 
therefore prefer in the following the form as given in (14). In [7] the difference 
between the two results was interpreted as an indicator for the magnitude of the 
"magnetic terms", that is, the contribution proportional to j. In the next section 
we will estimate the contribution of the transverse photons, making use of the long- 
wavelength limit. It remains to be seen, however, which one of the two possible schemes 
(longitudinal/transverse decomposition versus scalar/magnetic terms) will be better 
suited for explicit calculations. 

In order to determine total cross sections we would like to simplify the summation 
over all possible states. The expression (14) depends on w and A in two ways: 
Through the energy-conserving delta functions and through the expression for q 2 , 
where qf = A 2 /(/3 2 7 2 ) + uj 2 /(/3 2 j 2 ) + (2uA)/(/3 2 7) depends on both u and A. 
Replacing lo and A in qi by some average values u>o and Ao, we can perform the 
closure over all final states to get 

V^^Sinc,A(fc)S inc , n (s) (18) 
with q 2 , s 2 , and k 2 now the ones using ujq and Ao and 

Si„c,A(fc) = ^T|^0,A(fc)| 2 ( 19 ) 
X 

Si„c,n(s) =]T|F x , ji(s)| 2 . (20) 

X' 

In section 6 the dependence of the cross section on the choice of A and uj will be 
studied by varying both parameters over a reasonable range. For the atomic scattering 
function Sine, A, this study requires analyzing the contributions of individual electron 
excitations shell by shell, since the binding energies vary from some eV to several keV 
for the different shells. 



3. Contribution of transverse photons 



In section 2 we have only calculated the effect coming from the longitudinal photons 
(that is, coming from the charge operator). Already in [11, 7] the effect of the so-called 
"magnetic interaction" , that is, the effect of the current operator was estimated to be 
of the order of 1%. As the DIRAC experiment requires an accuracy of 1%, these 
contributions need to be considered more carefully. The part of the current operator 
j in the direction of q is already included in the above calculation. Therefore we 
need to study only the contribution coming from the transverse part of the current 
operator, that is, the effect coming from the transverse photons. Here we estimate 
them with the help of the long-wavelength approximation. Following [13], see also 
[14], W\ and W2 can be expressed in terms of the Coulomb and transverse electric 
(magnetic) matrix elements: 

Wi = 2tt (|T c | 2 + |T m | 2 ) (21) 

W 2 = ^2tt ^2|M c | 2 - (|T C | 2 + |T m | 2 )) , (22) 

where we use the usual definitions for M c and T°< m : W m = Att\M c \ 2 and W xx ' = 
<5 AV 27r (|T°| 2 + |T m | 2 ) with A, A' denoting the two transverse directions. The T m can 
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be safely neglected in our case. The long- wavelength limit relates T e to M c . In the 
multipole expansions 

T ° = E E t /m, ^ C = EE m ?m, (23) 

J>1 M J>0 M 

we have for a given multipole [8, 15] 

T C JM « j (^) 7 Mj C m, for J > 1. (24) 

In the case of the pionium only odd multipoles contribute and there is thus no M§ Q 
term. As higher multipoles are strongly suppressed we can safely set the factor 
( J + 1)/J to its maximum value at 2 (i.e. J = 1) for the purpose of an analytical 
estimate for the upper limit of \T e \. However, since our computer program already 
contains a multipole expansion, it is of course straightforward to perform a more 
refined numerical calculation. For the atom, on the other hand, Mq ^ 0; it is in fact 
the dominant contribution in the elastic process. Although generalizing the relation 
(24) with the maximum factor V2 for all multipole orders still provides an upper limit 
for T c , it may be a less useful over-estimate for the atom. The relation between |M C | 2 
and |T e | 2 is then approximately 

\T%q)\ 2 ^2^\M c (q)\ 2 (25) 

- (26) 

One sees that the transverse photons contribute to both W\ and W 2 . Denoting their 
contribution by Wj and W 2 , they can be expressed with the help of (21), (22), and 
(26) as 

W T = ^ W 0j (27) 

Wl = (28) 

First we look at the elastic case (for the atom). In this case we have A = 0, 
Wi^a — and W r f° A = |Foo,a(<z)| 2 . The product of the electromagnetic tensors is then 

W^W^n = -W^ n W 2 , A + l 2 WlnW 2 .A (29) 

= ^ (-l+7 2 ^)< f °nlWg)| 2 . (30) 

Again summing over all excited states of the pionium and neglecting the dependence 
on lo in q 2 and s 2 we get the total cross section 

4a 2 u? 



''-I 



i2 «^^ ("^4 + ?(fp)J l f »».A(9)l 2 - (3D 



We see that this cross section differs from (15) by the replacement 
1 u? 



s 4 s 2 y 7 2 g 4 + s 2(_ g 2) J ( 32 ) 

The estimate for the reduction of these terms compared to the charge contribution 
can be seen easily in this equation, if one assumes that the dominant contributions 
come from the range of q 2 (and therefore also of s 2 ) of the order of k^, where fcn 
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denotes the Bohr-momentum of the pionium, fc n = V a Bohr,n ~ 136.566/aeohr- This 
momentum is of the order to /a. Therefore the factor uj 2 / s 2 would give a reduction of 
the order a 2 , making this contribution completely negligible. However, a discussion 
of the relevant range of q 2 to be given in section 6 will show that this estimate is too 
crude. 

In the inelastic case (on the atom side) we can approximately set Ao~0 again, 
as the atomic binding energy is small compared to the other energies. Then the 
contribution from the transverse current of the atom will be suppressed even more 
than for the pionium, since the ratio A/k is even smaller than oj/s. The only difference 
compared to (31) is then the replacement of |F o,a(q , )| 2 by Si nc ,A(k). We get 

tLi = / d V^r^ ("^ + s 2 (-q 2 ) ) Sinc > n(s) 'WW' (33) 
We will discuss the contribution to the cross section in section 6. 

4. Dirac-Hartree-Fock-Slater model 

We now turn to the question how to evaluate the form factors and scattering functions 
derived in the general formalism of the preceding sections. For our specific application 
of pionium scattering off atomic targets, we use the pionium form factors as described 
in [7]; here we shall only discuss the calculation of the atomic form factors. 

The purpose of this explicit calculation of atomic form factors and especially 
atomic scattering functions is twofold: Existing tables [16] give no indication about 
the contributions from individual atomic shells, an information that is crucially needed 
to determine the appropriate average excitation energy in the closure approximation. 
Furthermore, since target atoms as heavy as Pt (Z = 78) will be employed in 
the experiment DIRAC, the atomic structure is more appropriately treated with 
relativistic orbitals, at variance with the non-relativistic calculations in [16]. 

4-1. Atomic ground state elastic form factors 

Within the framework of (Dirac-)Hartree-Fock-Slater theory, the atomic ground state 
wavefunction entering the expression for the form factor, 

z 

F on (k) = (* |^cxp(ifc.r,)|vI/ ), (34) 

is given by a single Slater determinant constructed from products of independent 
particle orbitals, 

*° = "^E si S n WV)( r i)'" $ p(z)( r 2)- (35) 

v ■ p 

Here Z denotes the nuclear charge (and the number of electrons), while p denotes the 
permutations of orbital indices, and $j signify single particle orbitals. Each of the Z 
exponential terms exp(ifc • rj) in (34) acts as a one-particle operator. Orthogonality 
of the orbitals effectively cancels the summation over permutations in the bra and ket 
vectors leading to 

z 

Fov(k) = 5>,(r,)| exp(ife • r,)|$,(r,)). (36) 
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So far we have not considered any angular momentum coupling of the independent- 
particle orbitals, that is, our ty is determined by a set of quantum numbers 
(nj,lj,m,j), j = 1 . . . Z, without coupling to a total L (or J) and M. Furthermore we 
assumed the ground state wavefunction in the bra and ket symbols to be identical. 
This need not actually be the case, as different z-projections M of the total angular 
momentum J of the ground state cannot be distinguished (without external fields). 
The angular momentum coupled ground state wavefunction would then be obtained 
by summing and averaging over M and M' in the bra and ket vectors, respectively. 
Likewise on the level of independent particle labels nij, the orthogonality argument 
applies strictly only to all orbitals except j, that is, we should distinguish between 
nij and m'j for the bra and for the ket vector, respectively. Since the one-particle 
operator does not affect the spin part of the wavefunction, we should also insert an 
orthogonality factor for the spin orbitals in bra and ket, Xj, x'j- Of course the principal 
and azimuthal quantum number nj and lj coincide in the bra and in the ket symbol. 
Expanding the exponential in spherical harmonics we find immediately 
z 

F 00 (fe) = ]T(2^ + f)<5 Xj , x ^i^y A * M (fc)V2ATT 



i=i a,^i 



A 



l 3 l 3 





m,j —nij t 1 
with the radial form factor defined by 

2 



H?j(k)= drr z R nih {r)j x {kr)R n]h {r), (38) 
Jo 

where the R n i{r) denote radial wavefunctions for the orbitals, and j\(kr) is a spherical 
Bessel function. 

In a next step averaging |Foo(fc)| 2 over all directions q and using the orthogonality 
relation of the spherical harmonics yields 



|F 00 (fc)| 2 :=^ J dfc|F 00 (fc)| 5 



z 



= £(2A+i) m-ir^+i)<5 : ' ' ' '- ; 



lj lj A 
x 3 ,r 3 V 



A,ju ( 3 = 1 

' 7^(fc)' - (*M 



lj lj A 



nij —nij n 

Obviously all electrons contribute coherently to the form factor, as expected. Due to 
the first 3j-symbol only even multipoles contribute to the sum. 

In the LS-coupling scheme, the atomic ground state is characterized by a specific 
total L and a (possibly averaged) total M. Instead of coupling the individual electrons' 
angular momenta to total L and then averaging over M, we average directly over 
individual m,j,m'j for the orbitals occupied in accordance with the Pauli principle. 
This amounts to neglecting energy differences between fine structure levels. Hund's 
rules (see e.g. [17]) state that sub-shells are to be filled by adding as many electrons 
with different nij and the same spin projection as possible. The critical multiplicity 
is thus that for a half filled sub-shell, (21 + 1). For half filled or completely filled sub- 
shells, Hund's rules imply that both m,j and m'j run over all possible values from — lj to 
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Consequently the angular part for these spherical sub-shells reduces to selecting 
only monopole contributions. A completely filled sub-shell (rij,lj) thus contributes 

2(2i j + i)n° jj (k) 

to the full form factor. (The leading factor 2 indicates the spin multiplicity.) 

For sub-shells that are neither completely nor half filled the averaging procedure 
yields different multiplicities depending on the multipole order. For the dominating 
monopole contribution, this factor is simply given by the occupation number of the 
sub-shell, whereas for higher multipoles this factor is proportional to the product 
of occupation number and the number of holes in a sub-shell with identical spin 
projections. For a given number of electrons in an open sub-shell we averaged the 
m-dependent part in (39) over all possible distributions of rrij and m'j values. 

Note that the coherent (elastic) form factor F 00 (k) derived above describes only 
the effect of scattering off the atomic electrons. The complete elastic form factor for 
the atom reads 

^Atom(fc) =Z-F 0Q (k), 

assuming a point-like nucleus. 

4-2. Atomic inelastic scattering functions 

Besides the elastic form factor F m (k) treated in the preceding section, we also need 
to consider the contributions due to excitations of the atomic electron cloud. Nuclear 
excitations will not be considered here because the much larger excitation energy 
required (typically on the order of MeV) exceeds the energy range relevant for our 
application to pionium-atom scattering. We demonstrated in [7] that deviations from 
a point-like nucleus are negligible for the electromagnetic processes considered here. 
Thus the nucleus' internal structure with its excited states is equally irrelevant as the 
experiment DIRAC cannot probe this structure. In analogy to the elastic form factor, 
a transition form factor is written in the form 

z 

F XQ {k) = (*x|5>xp(ifc-r,)|* ), (41) 

corresponding to scattering with excitation of the atomic electrons from the ground 
state to some excited state X. A similar expression was studied in [18] in the context 
of the equivalent photon approximation. The total inelastic scattering function is 
defined as the incoherent sum over all states X other than the ground state 

= ^|^o(fc)| 2 -|F 00 (fc)| 2 

all X 

Z 

= Z + Y^{M cx P (ife • [r, - ri])|tt ) - l^oo(fc)| 2 . (42) 

Here we have added the ground state in order to exploit the completeness of the set 
of states {X}: Expanding the squared modulus of Fxo introduced a second (primed) 
set of variables r'j which has been removed again by virtue of the completeness of the 
set of states {X}. Furthermore we evaluated the sum over the diagonal terms i = j 



(40) 
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separately, obtaining the term Z (since exp(ife • [rj — r,]) = 1 in this case). Using the 
same Slater determinant wavefunctions for as in the preceding section, expanding 
the double sum corresponding to |F 00 (fc)| 2 , and combining terms with the last sum, 
we find (in terms of the single-electron orbitals) 

z z 

S inc (k) = Z-J2 ^ |(^| exp(ife • r)|$,)| 2 . (43) 

i=l 3 = 1 

The new terms required in the determination of 5*i nc are the matrix elements 
($,|cxp(ifc • r)\^j) = W-irV 4 ^ + 1)(2 ^ + l)E y A, M (fc)i A V2ATT x 

( o s o ) ( x x i ) <«> 

some of which (namely, those with i = j) have already been used in calculating 
F o(fe). As before, XiiXj denote the spin projections. Averaging over the directions k 
we obtain immediately 




Z Z 2 

= ^-EE^x 3 (^ + l)(2^ + l)E( 2A + 1 )( o o o) X 

i=lj=l A V ' 

( h h X ) 2 [n*(k)] 2 . (45) 

Considering again the non- averaged incoherent scattering function in the special case 
of completely filled sub-shell k, lj, the summations over m, and mj remove the angular 
dependence on k. The filled sub-shells thus contribute 

2(2Z i + l)(2Z J . + l)]T(2A+l)( |j \ o) 2 [^.(fc)] 2 

to the incoherent form factor, with a factor 2 for the spin multiplicities in both sub- 
shells (rather than a factor 4, because the cross terms between sub-shells with opposite 
spin projections drop out due to the orthogonality of the spin orbitals). 

Except for Z = 58 (Ce, not of interest to us) all atoms with Z < 90 have in 
their ground state only one sub-shell that is neither half nor completely filled [19]. (In 
atoms with two open sub-shells, one of them has I = 0.) Thus we need not consider the 
m-averaging procedure for cases where rm and rrij come from two different partially 
filled sub-shells with k > and lj > 0. In the cases of our interest averaging over 
rm and mj is then straightforward. The m-averaged contribution from two different 
sub-shells reads 

\, iV ^(2X + l)^ \ J) 2 [7^.(fc)] 2 

where vi and i/j denote the occupation numbers of the two sub-shells and the factor 
1/2 stems from the orthogonality of spin orbitals. If i and j both refer to the same 
sub-shell the multiplicity factor is slightly more complicated, depending on whether 
the sub-shell is less than, or more than, half filled. For this limited number of cases, 
we again determined the m-dependent part by suitably averaging over all possible 
distributions of m values in a partially filled sub-shell. 
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5. Some other models 



For comparison we briefly discuss in this section some simple alternative models for 
evaluating the coherent and incoherent form factors, Foa(k) and Si nc (k), respectively, 
for application in complex atomic scattering. Specifically, we will use analytical 
screening models [20, 21] to derive the elastic form factors. The inelastic scattering 
functions can then be obtained from the elastic form factors either in the no-correlation 
limit, or from an argument due to Heiscnberg and based on the Thomas-Fermi model. 

The simplest possible model to describe the effect of incoherent scattering off the 
atom's electrons would merely divide the cross section for coherent scattering (scaling 
with Z 2 ) by Z, on the grounds that everything remains the same except that each 
electron contributes incoherently to the cross section (complete "anti-screening"). We 
found that this approach underestimates the incoherent scattering cross section by as 
much as 50%. For typical targets like Ti or Ni this implies an error of roughly 2% in 
the target-inclusive cross section, clearly beyond the required limit of 1%. 



5.1. Elastic form factors 

In order to simplify the atomic structure calculation, one might use the Thomas- Fermi 
model to replace the density p(r) occurring in 

F a{k) = J d 3 rp(r) cxp(ife • r). (46) 

Expressing the potential due to the charge distribution of the electrons in the form 

V(r) = -f xW, (47) 

the corresponding charge distribution is given by the second derivative of the screening 
function: 

P(r) = ^- r X"(r). (48) 

Here and in the following the prime denotes differentiation with respect to r. For a 
spherical charge distribution the coherent form factor reduces to the monopole term. 
Using Molicrc's [21] parameterization for x, 

3 

x(r)=^B l exp(-Ar/6); (49) 
i=i 

B 1 = 0.1; B 2 = 0.55; B 3 = 0.35; (50) 
j3 1 = 6.0; 02 = 1.2; fa = 0.3; (51) 

with b = aBohr(97r 2 /128) 1 / 3 Z -1 / 3 , the coherent form factor reads 

^<*>-*2>ITI5W (52) 

for the electronic part, and again -FAtom(^) = Z — Foo(k). The same analytical 
form with different parameters B t and Pi determined by fitting expectation values 
of powers of r to exact Dirac-Hartree-Fock-Slater results is used in [20] where the 
fitting parameters for all Z < 92 may be found as well. 
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5.2. Inelastic form factors: No- correlation limit 

Inserting the elastic form factor of the previous subsection into (42) we are left with 
the evaluation of 

z 

exp ( ifc ' fa ~ r *H*o) = Z{Z - 1) / d 3 r dV JV 2 (r,r') exp(ifc • [r - r'])(53) 

i=i j^i 

where we have already integrated over all variables not pertaining to the orbitals i 
and j. Here the function N 2 (r, r') describes the probability of finding any two of the 
properly anti-symmetrized electrons at positions r and r' . The no- correlation limit 
now replaces this two-particle density by the product of the single-particle probabilities 
p(r)/Z and p(r')/Z. The integral in the last expression then reduces to \F 00 (k)/Z\ 2 , 
i.e., the square of the elastic form factor normalized per electron. As there are Z(Z—1) 
such terms in the double sum over i,j, we finally find 

S inc (k) = Z- \F 00 (k)\ 2 /Z. (54) 

Using the elastic form factor of the previous subsection, this result provides a simple 
expression for S[ nc . However, we have made here the crucial assumption that there is 
no correlation between the electrons in different orbitals. In a single-particle picture 
this amounts to assuming that all single-particle states are available for all electrons 
simultaneously. Pauli blocking, i.e., the fact that due to the Pauli exclusion principle 
some states X cannot be excited for a given electron because they are occupied by 
other electrons, is then disregarded completely. Since in this limit the (incoherent) 
summation over X also includes the Pauli blocked states, the expression (54) clearly 
overestimates the correct scattering function. 



5.3. Thomas-Fermi model for incoherent scattering 

Following Hciscnbcrg [22] we can find a simple expression for the incoherent atomic 
form factor, going beyond the no-correlation limit but remaining in the spirit of the 
Thomas- Fermi model (see also [23] ) . Expanding the squared modulus in the incoherent 
sum (43) we write 

z 

S inc (k) = Z- / d 3 rd 3 r' exp(ifc- (r-r'))|^ $ j( r ) $ j( r ')| 2 - (55) 

In the Thomas-Fermi model, the density is related to the volume of a sphere in 
momentum space, 
z 



E I'W 2 = T^W / d " K = i[2|^(r)|] 3 / 2 . (56) 

Heisenberg generalized this expression to obtain the two-particle density as an integral 
in momentum space (with n = (r + r')/2) 

]T $*(r)$,(r') = / d 3 « exp(i« • (r - r')) (57) 

which reduces to the Thomas-Fermi expression for r = r'. (55) now contains 
additional six-fold integration over d 3 K and d 3 K ; . Performing these integrations as well 
as the integration over (r — r'), we find (following Heisenberg's argument about the 
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common volume of two intersecting spheres in momentum space, and paying attention 
to the spin multiplicity) 

4 f r ° 

S htc (k) = Z-—J dr 1 rl(^m^-k/2)\^/m^) + k/A) 1 (58) 

where the integration over n is restricted to the region of coordinate space with 
fcr(^i) > k/2 since otherwise the two spheres in momentum space do not overlap. 
Using V(r) = (Z/r)x(r) the incoherent form factor turns into 

4 f r ° 2 f r ° 1 

S inc (k) = Z-— &r^[2Z X {r)f /2 + -Zk dr r X (r) - —(kr ) 3 (59) 

where the upper limit of integration r must satisfy 

-XM = \k 2 . (60) 

Within the frame work of the Thomas-Fermi model, the screening function 
satisfies the differential equation 

f^X(r) = ^[(2Z/r) X (r)f 2 : x (0) = 1, (61) 

thus enabling us to replace the first integral in (59). The second derivative of \ can 
then be removed by integrating by parts, yielding 

2 f r ° 1 1 

S inc (k) = -Zr x'(r ) + -Zk dr r X (r) + -r k 2 - ^(krof. (62) 

TT J Q 8 60TT 

From (59) and (62) we note that this simple model does not reproduce the correct 
limit as k — > 0: In this limit r — > oo, and the integral in (62) assumes a constant 
non-zero value. Thus S lnc grows linearly with k. By contrast the expression in the no- 
correlation limit (containing F 00 (k)) grows with k 2 , as does the Hartree-Fock-Slater 
result derived in the preceding section (since the term linear in k drops out of the 
expansion in (42) due to the symmetry under interchange i <-> j). 



6. Numerical method and results 



In our application to scattering of pionium on atomic targets, we use hydrogenic 
wavefunctions for the pionium system 7r + 7r~. The corresponding form factors F o,n 
and S'inc.n are evaluated analytically as described in [7]. 

Section 4 provides expressions for the evaluation of coherent atomic form factors 
and incoherent scattering functions in the framework of the Dirac-Hartree-Fock-Slater 
formalism. In our calculations we start from a simple analytical charge distribution for 
the electrons as given e.g. in [20] or similarly in [21]. Taking the electronic structure 
of the elements from [19] we solve either the Schrodinger or the Dirac radial equation 
for each occupied orbital, treating exchange effects by Latter's approximation (see 
e.g. [20]). The resulting charge density is then used to obtain improved radial orbitals, 
iterating the process until self-consistency is reached. Even for heavy elements with 
electrons in some twenty different orbitals and requiring several ten iterations, this 
calculation is readily performed with the help of the program package RADIAL [24] . 

Using these orbitals we evaluate the radial integrals in (39) and (45) on a 
reasonably dense mesh of k values with the help of an integration routine developed 
for integrals containing spherical Bessel functions and powers [25]. To this end, the 
numerical solutions for the orbitals obtained on a grid of r values are replaced by 
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Figure 2. Various contributions to the integrand for the cross section vs. qj_, 
on a log-log scale. To compensate for the logarithmic gx-axis, the integrand is 
represented as da/dlnqj_. The individual curves are labelled in the figure. The 
cross sections correspond to target-elastic (solid lines) and target-inelastic (dashed 
lines) pionium scattering off Ni (Z = 28) at a projectile energy E = 5 GeV, 
summed over all final states of the pionium. The arrows at the upper edge indicate 
the relevant momentum scales. Sec text for details. 



piecewise splines. The angular parts for the partially filled sub-shells are determined 
by averaging over all distributions of magnetic quantum numbers m in accordance 
with the Pauli principle, as described in section 4. 

Let us first investigate the range of q± relevant for the cross section as given by (9) 
or (18). Figure 2 demonstrates the interplay of different momentum scales associated 
with the photon, with the atom, and with the pionium, respectively. The figure shows 
the integrand from (9) in singly differential form: 



(Photon) x (Atom) x (Pionium), 



dcr Aa 2 

d^I ^ 

together with its decomposition into photon, atom and pionium parts (cf. (8)): 



(Photon) = — 



1 



(Atom) = W 2 ,a 



I^ocaWI 2 

(Pionium) = W 2 ,u = (q A /s 4 )S incM (s). 



for coherent scattering, 
for incoherent scattering, 



(63) 



(64) 



(65) 
(66) 



The solid lines refer to the total cross section (i.e., summed over all pionium final 
states) for target-elastic scattering off Ni (Z = 28) for pionium in its ground state. 
The projectile energy is 5 GeV. The dashed lines correspond to the same setting for 
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the target-inelastic process. Here the integrations over u and A in (9) have been 
replaced by setting an average excitation energy for the pionium at ujq — 1.858 keV 
(ground state binding energy), and for the atom we set A = (target-elastic) and 
A =100 cV (target- inelastic), respectively. 

The arrows indicate the relevant momentum scales: qi for the photon, /ctf = 
-Z^/aBohr for the atom, and ku for the pionium. The arrow on the left under the label 
u qi" corresponds to target-elastic scattering, while the one on the right corresponds 
to incoherent scattering (with non-zero Ao and thus with a larger qi). For qj_ <C qi, 
the photon momentum is essentially given by the constant qi . When q± ~> qi, on the 
other hand, this part displays a l/q]_ behavior. The atomic part shows an increase 
between qi and several inverse Bohr radii (indicated by /ctf)- As expected, |-Fbo,A(fc)| 2 
grows roughly with q\, whereas Si nCy A(k) grows only with q\. At q± w 5&tf the 
atomic part reaches its asymptotic value (Z 2 or Z, respectively). In this regime the 
pionium part only just starts to contribute appreciably. It grows quadratically with 
q± to saturate at a few multiples of the pionium scale fen- The product of the three 
factors clearly demonstrates that the main contributions to the cross sections come 
from the range of q± between fcxF an d fen- 

In figure 2 we set luo — i^bind.n for the pionium, as well as a non-zero value 
for A in the case of incoherent scattering. The specific choice for these average 
excitation energies is guided by the following observations. For coherent (target- 
elastic) scattering A = 0, and any ambiguity in calculating total cross sections from 
(18) is limited to the choice of an appropriate wq. From a comparison [7] of total cross 
sections for coherent scattering in the closure approximation with the "exact" total 
cross sections obtained by accumulating partial cross sections for bound-bound and 
bound- free pionium transitions (summation/integration over all final states), we note 
that bound-bound transitions (excitation and de-excitation) account for the major 
part of the total cross section. Typically, breakup (ionization) accounts for some 
30%-40% of the total cross section in the ground state, decreasing roughly by a factor 
n for pionium in the initial state (n, I). Furthermore, most of the breakup cross 
section from a given initial state comes from the range of continuum energies from 
to about i?bind,n above the continuum threshold. Therefore the average energy 
difference between initial and all final states — weighted by their contribution to the 
total cross section — is of the order of the binding energy -Ebind.n- 

Figure 3a) shows the total cross section for elastic scattering in the closure 
approximation as a function of loq (in units of -Ebind.n)- The total cross sections 
have been normalized to their values at loq = 0. As can be seen from the figure, only 
the ground state total cross section varies appreciably over a reasonable range of Uo 
values. We also find [7] that the closure cross sections at ui = -Ebind.n coincide with 
the converged accumulated partial cross sections. We therefore set ui a = i?bind,n in 
the following. 

In frame b) of the same figure we show the dependence on the average atomic 
excitation energy Ao for the target-inelastic scattering process. Here the pionium 
energy has been fixed at ujq w 1.858 keV. Since we calculate S- mc {k) using individual 
orbitals, we can easily determine the contributions resolved with respect to atomic 
shells, or even per individual electron as in figure 3b). The tabulated values for 
the atomic incoherent scattering functions [16] cannot serve to study the closure 
approximation's sensitivity to the atomic excitation energy. When varying Ao from 
to several keV, it is important to distinguish whether this excitation energy is 
transferred to an electron in the K shell or to one of the outer electrons. In the latter 
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Figure 3. a) Closure cross sections (normalized to loq = 0) versus u>o (in units of 
^bind.n) for coherent scattering of pionium in various initial states as indicated in 
the figure. Target material is Ni, projectile energy is 5 GeV. Solid lines: initial s 
states; dashed lines: initial p states of the pionium. b) Contributions of individual 
target electrons to the closure cross section for incoherent scattering of pionium in 
its ground state, plotted versus average atomic excitation energy Ao (normalized 
to the average binding energy of the respective atomic shells, ranging from some 
10 eV for the ./V-shell to ~ 5 keV for the i-sT-electrons) . Also shown is the average 
contribution of all shells (total incoherent cross section divided by Z). Target 
material is Ti, projectile energy is 5 GeV. 



case, the atomic electron is excited into a high energy continuum state. This process 
contributes very little to the incoherent cross section which — like the pionium's cross 
section — is dominated by excitation rather than by ionization. The average excitation 
energy A for each curve in figure 3b) has been normalized to the binding energy of the 
individual shells (averaged over sub-shells). We note that the individual contributions 
are roughly proportional to the principal quantum number of each electron. Combined 
with the fact that the outer shells typically accommodate many more electrons than 
the inner shells, we find that the target-inelastic (incoherent) scattering process is 
clearly dominated by the loosely bound outer electrons. At the same time figure 3b) 
also demonstrates that the incoherent scattering cross section is almost independent 
of Ao. Thus we may safely set Ao = in our calculation. 

As a verification and further illustration of these findings, figure 4 displays 
dai nc /dlnq± for ground state pionium scattering incoherently off Ti (Z = 22) at 
5 GeV projectile energy. The solid line corresponds to the integrand for incoherent 
scattering (18) calculated using (45). In these calculations we set ujo = 1.858 keV, 
the binding energy of the pionium, and A = 0. As can be seen from the dashed 
line in the figure, the simplest approximation consisting of scaling the target-elastic 
(coherent) cross section by 1/Z clearly underestimates the correct result. At the same 
time the dash-dotted (chain) curve shows that the approximation using the coherent 
form factor (39) in the no-correlation limit obviously overestimates the correct result 
by an even larger amount. In more detail we find that the ratio (iJcoh/Z) : o-- mc for 
pionium initially in the ground state amounts to 0.49 (!) for a Li target, 0.62 (Al), 
0.66 (Ti), 0.73 (Ni), 0.72 (Mo), and 0.74 (Pt). For pionium initially in 2s or 3s 
states, these ratios range from 0.65 (Li) to 0.80 (Pt) for these targets relevant in the 
experiment. For the lighter targets, these huge uncertainties of the incoherent part 
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Figure 4. Integrand for incoherent (target-inelastic) cross section of ground 
state pionium scattering off Ti at energy 5 GeV, versus q± (on a log-scale). 
The areas under the curves yield the respective cross sections. Calculations with 
uiq = 1.858 keV and Arj = 0. The contributions from the individual shells (dotted 
lines) are labelled in the figure. 



result in considerable errors in the inclusive cross sections: The ratio 
a coh (l + l/Z) 

&coh ^inc 

amounts to 0.795 (Li), 0.958 (Al), 0.978 (Ti), 0.987 (Ni), 0.991 (Mo), and 0.996 (Pt). 
Thus only for the heavy targets the required accuracy of 1% can be attained with such 
a crude approximation for the incoherent part. 

Further in figure 4, the dotted lines show the contributions to the incoherent 
scattering cross section resolved according to atomic shells. A direct comparison of the 
areas enclosed between the individual lines and the abscissa demonstrates that the ten 
M-electrons dominate the cross sections, followed by the eight L-electrons. Also, the 
two iV-electrons contribute considerably more strongly than the if-electrons, whose 
influence is limited to rather large q± as expected. 

Figures 2 and 4 show clearly that the principal contributions to the cross sections 
come from the region /ctf < Q± < ku- Thus the long- wavelength limit applies for the 
pionium (but not for the atom). Furthermore we noted that the relevant excitation 
energy A is small and may safely be set to zero, and thus Wi.a = 0. Under these 
circumstances the cross section due to the transverse part of the current operator is 
given by (31) and (33). In figure 5 we show the ratio between the cross sections for the 
transverse electric part (T°) and the one for the Coulomb part (M c ), as a function 
of the average pionium excitation energy ujq. Note that the ratio is given in per cent. 
The solid lines correspond to total cross sections for coherent scattering of pionium 
in various initial s-states as indicated in the figure. The dashed lines correspond to 
the target-inclusive process (with Ao = 0). At the typical average pionium excitation 
energy uj = -EWd the transverse cross section cr T for coherent scattering amounts to 
0.4% of the Coulomb part in the ground state, decreasing rapidly for initially excited 
states of the pionium. 
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Figure 5. Ratio (in per cent) between transverse electric and Coulomb 
contributions to the total cross section of pionium scattering off Ni at 5 GeV, 
as a function of luq. Solid line: target-elastic process; dashed line: target- inclusive 
process. Pionium initial states (with Zj = 0) are indicated in the figure. 

Finally figure 6 compares various models for -Foo,A and Si nc ,A over the range of 
k values relevant for our calculations of cross sections in the context of experiment 
DIRAC. The solid lines correspond to our calculations using (39) for the coherent 
form factor and (45) for the incoherent scattering function, respectively. We use 
relativistic Dirac orbitals in both cases. The squares represent tabulated values from 
two compilations of state-of-the-art Hartree-Fock calculations with various corrections 
(configuration interaction, relativistic effects, a.s.o.). Note that while the tabulated 
results for Fqq correspond to relativistic calculations [26], the tables for Si nc in [16] 
contain only non-relativistic results. The log-log representation in the left frame serves 
merely to display the asymptotic behavior of the simple models discussed in section 5. 
Analyzing the dashed line corresponding to (52) we note that the incorrect asymptotic 
fall-off for Fqo at large k poses no difficulties. In the right panel we see that between 
5 < k < 10 a.u. the Thomas-Fermi-Moliere model misses features of the atomic 
shell structure, but the resulting deviation from (39) is insignificant. Much more 
problematic are the crude approximations for 5*i nc . While the no-correlation limit 
(54), using (52) and shown with the dotted line, increases with k 2 at small k (as it 
should), it dramatically overestimates the scattering function in the most relevant 
range between 0.1 and 100 a.u. On the other hand, the Thomas- Fermi model for 
incoherent scattering as developed by Heisenberg [22], with modifications [23, App. 
B] for the use of Moliere's approximation, is quite successful at k > 5 a.u., but it fails 
completely at smaller k. 

7. Conclusions 

We have reviewed the formalism for incoherent atomic scattering. The basic 
expressions for the atomic form factors and scattering functions have then been 
evaluated in the framework of Dirac-Hartree-Fock-Slater theory, i.e., using numerically 
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Figure 6. Electronic part Foo of the coherent atomic form factor and incoherent 
scattering function Si nc for Ti (Z = 22). The asymptotic behavior is more easily 
seen from the log-log diagram on the left. The range of relevance for the cross 
section calculations is 0.1 < k < 100 a.u. For an explanation of the different 
models, see text. 

determined electron orbitals. For comparison, both the form factor for coherent 
scattering, as well as the incoherent scattering function have been derived in simple 
analytical models based on the Thomas-Fermi model of the atom. Applying these 
different descriptions we performed detailed numerical studies in the context of 
pionium scattering incoherently off the electrons of various target atoms. Due to 
the much larger reduced mass of the pionium system (^ n = to tt/2 ~ f36.566m e ), the 
length and momentum scales of the pionium and the target atom are very different. 

An investigation of the relevant momentum transfer revealed that the cross 
sections are dominated by the contributions from the region between the Thomas- 
Fermi momentum /ctf = •^ 1 ^ 3 /aBohr for the atom, and the momentum scale of 
the pionium at fcn = Mn/oBohr- Under these circumstances the simple models for 
incoherent scattering discussed in section 5 prove not sufficiently accurate for our 
application in the context of pionium-atom scattering with a required accuracy of 1% 
or better. Whereas the analytical models discussed in section 5 provide sufficiently 
accurate coherent form factors, the incoherent contribution really requires the more 
accurate treatment developed in section 4, despite its lesser importance (as compared 
to coherent scattering). Only an explicit DHFS calculation can provide satisfactory 
scattering functions. Since the target materials employed in experiment DIRAC will 
be as heavy as Pt, our calculation with relativistic orbitals is clearly more appropriate 
than the non-relativistic results tabulated in [16]. 

From our detailed discussion of the ^-dependence of the integrand for the cross 
section we also conclude that the loosely bound outer shell electrons dominate the 
target-inelastic cross sections. Their contribution to the integrand dcr/dg_i_ covers a 
much larger range of q± values than the one corresponding to inner shell electrons. 
Following this argument further, free electrons would show a behavior rather similar 
to the one of the outer shells, their contribution to the cross section would stretch 
even further down to smaller qj_. However, the cross section hardly depends on this 
modification at very small q± . Our calculation therefore applies equally well to quasi- 
free electrons and to electrons in the conduction band. Thus solid state effects and 
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chemistry need not be considered explicitly as they prove irrelevant for calculations 
pertaining to experiment DIRAC. The same conclusion had been inferred in our earlier 
study [7] based on an analysis of the relevant impact parameter range. Recalling that 
target-inelastic scattering constitutes merely a correction of order 1/Z of the atomic 
part, variations on the order of a few per cent in the incoherent scattering cross sections 
are insignificant. 

These findings are confirmed in our analysis of the dependence of the cross 
sections on the average excitation energies for pionium and atom. Our calculation 
of incoherent scattering contributions resolved according to individual target electron 
shells demonstrates that the average excitation energy for the atom may safely be 
set to zero, Ao = 0. This information which is crucial for the closure approximation 
cannot be extracted from the scattering functions in [16]. For the pionium, on the 
other hand, a non-vanishing excitation energy in the amount of the binding energy is 
needed when calculating total cross sections in the closure approximation in order to 
get agreement between the closure result and the result obtained by explicitly summing 
the partial cross sections over all final states. 

Earlier investigations on pionium-atom interaction [11, 7] invoked properties of 
the hydrogen-like pionium system to obtain crude estimates for the magnitude of 
the magnetic terms so far neglected in this interaction. From the non-relativistic 
wavefunctions for the pionium, its internal velocity is of the order of v/c « a/2. 
Thus magnetic terms are believed to be small of this order. We established a much 
better justified estimate for the magnetic terms (actually for the transverse part of 
the current) in the long-wavelength limit. Our investigation showed that this limit 
applies very well for the pionium, whereas it does not apply for the atom. However, 
the transverse current does not contribute in the elastic case on the atom side, and 
for the inelastic case it is suppressed even more than the corresponding term on the 
pionium side because the relevant atomic excitation energies (of the outer shells) are 
smaller than those of the pionium. 
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